Spin injection into a short DNA chain 



in 
o 
o 

(N 
o 

Q 

0\ 



o 



T3 

C 

O 
o 



> 

in 
o 

-I— > 



-a 
c 
o 



X 
S3 



X. F. Wang and Tapash Chakraborty 
Department of Physics and Astronomy, The University of Manitoba, Winnipeg, Canada, R3T 2N2 

Quantum spin transport through a short DNA chain connected to ferromagnetic electrodes has 
been investigated by the transfer matrix method. We describe the system by a tight-binding model 
where the parameters are extracted from the experimental data and realistic metal energy bands. 
For ferromagnetic iron electrodes, the magnetoresistance of a 30-basepair Poly(G)-Poly(C) DNA 
is found to be lower than 10% at a bias of < 4 V, but can reach up to 20% at a bias of 5 V. In 
the presence of the spin-flip mechanism, the magnetoresistance is significantly enhanced when the 
spin-flip coupling is weak but as the coupling becomes stronger the decreasing magnetoresistance 
develops an oscillatory behavior. 

PACS numbers: 87.14.Gg,72.20.Ee,72.25.-b 



In recent years, a remarkable progress in direct mea- 
surements of electron transport through DNA has gen- 
erated intense interest in DNA electronics 0, 0, y[ 0, 
HISIZI,EH[I3[HE0- DNA is found to have 
diverse electronic properties depending on its structure 
and the environment around it The clear semicon- 
ductor behavior observed in a short DNA chain of 30- 
basepair Poly(G)-Poly(C) has been explained by a tight- 
binding model 0i 0, [3- On the other hand, spin trans- 
port through nanostructures has, of late, been receiving 
considerable attention because of the possiblity of devel- 
oping spin-based electronic devices [l4|. Inspired by the 
broad interest in spin-injection into mesoscopic systems 



J 



h5|, we have investigated the quantum spin transport 
[l(j through a short DNA chain connected to ferromag- 
netic electrodes. We predict an enhancement and os- 
cillation of magnetoresistance in this system taking into 
account the realistic band structure of ferromagnetic Fe 
electrode and a spin-flip mechanism. 

We consider a p-type semiconductor DNA chain of N 
basepairs connected to a circuit via metal electrodes. The 
tight-binding Hamiltonian of the system is 



H — Hd + Hl + Hj 
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Here a is the creation operator of electron with spin 
a on site n {— 1, • ■ ■ , N) of the DNA chain, the left elec- 
trode (n < 0), and the right electrode (n > N + 1). Hd 
describes electrons (holes) of spin a in the DNA chain 
with the on-site energies —e d which is equal to the 
highest occupied molecular orbit (HOMO) energy of each 



base-pair, and the hopping parameters t" d between neigh- 
boring sites. The HOMO energy band is then determined 
by e a d and tZ. The Fermi energy Ep in the p-type DNA 
chain locates between the HOMO and lowest unoccu- 
pied molecular orbit (LUMO) bands and is closer to the 
HOMO band edge. Experimental results have indicated 
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FIG. 1: Energy band of the system in equilibrium. 



that the Fermi energy may vary from sample to sample 
2j. H x with X = L or R denotes the Hamiltonian of 
electrons in the left electrode (L) or the right electrode 
(R). In the tight-binding model, e'^ nX is the center of the 
energy band where the electrons are in the metals and 
^mX its band width. When the DNA chain contacts to 
the metal electrodes, exchange of electrons (holes) be- 
tween the DNA chain and the electrodes becomes pos- 
sible. In equilibrium, as illustrated in Fig. ^ the Fermi 
energies of the electrodes and of the DNA match and 
a tunnelling barrier forms between them. The contact 
property is described by the tunnelling parameter tdm- 
When a bias voltage drop is applied over the electrodes, 
distribution of the voltage drop or the potential profile 
along the non-equilibrium system depends on the DNA 
chain property and its contact with the metal electrodes. 
Since the free electron density in the metals is much 
higher than that in DNA, we assume that the band struc- 
ture of the metal electrode is not changed by the applied 
bias. When the contact between DNA and the electrodes 
is poor, the voltage drop concentrates on the contact. In 
case of a perfect contact however, the whole voltage drop 
should be applied mainly along the DNA chain. In this 
letter, we assume the voltage drop is on the contact since 
it is supported by the fit to the experimental result (see 
below) . 

In real world, a DNA chain is composed of two strands 
of bases with one phosphate-sugar backbone connected 
to each strand. The backbones can affect the on-site 
energy of electrons in the basepairs. Further, the envi- 



ronment around the DNA chain may also play a role in 
the property of the electrons. Here we use a reservoir of 
semi- infinite chain 0,0 03 

with a energy band of width 
47 and a coupling of strength ij to each DNA basepair 
site to mimic the effect of the backbones and the envi- 
ronment. As a result, the on-site energy of each site in 
the DNA basepair is modified by a self energy E£(i?) 
which is energy dependent and is expressed as, 



E 



with e r being the on-site energy of the semi-infinite reser- 
voir chain which we assume to be equal to the DNA on- 
site energy and S r = (E - e r )/2 — i[j 2 - (E - Erf/A} 1 / 2 
the self energy of any reservoir site which is obtained 
self-consistently. In what followSjWe have used the val- 
ues 77 = 0.1 eV and 7 = 5 eV [3|. To study the spin 
relaxation for a possible spin injection, we introduce the 
term H sp in the Hamiltonian to take into account the 
spin flip on-site and between neighboring sites described 
by the parameters t s £ and t s d respectively. The spin-flip 
along the DNA can be a result of spin-orbit interaction, 
magnetic impurity in the backbone, or magnetic environ- 
ment. 

In order to evaluate the transport properties of the sys- 
tem, we have employed the transfer matrix method 0,0]. 
For an open system, the secular equation of the system 
is expressed as a group of infinite number of equations of 
the form, 



The wave functions of sites n + 1 and n are related to 
those of sites n and n — 1 by a transfer matrix M, 
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and A n ,„ + i = (i„,„+i) 2 - (t s n n+1 ) 2 . Assuming plane 
wave functions for the electrons $„ = Y^ a (A a e lkLna + 



B° e - ikLna ) for 71 < and $„ = £ CT C a e ikLna for n > 
N + 1 in the left and right electrodes respectively, we can 
express the output wave amplitude C a in terms of the 
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FIG. 2: Energy dependence of the parameters e m (b) and 
t m (c) for ferromagnetic Fe. Solid curves and dotted curves 
correspond to spin-up and spin-down electrons. The resulting 
bulk DOS is also shown in (a). 



input wave amplitude A a and the transmission, 

|C ± | 2 sin(fe|a) 



T (E) 



|A+| 2 sin(fc+a) + |,4~| 2 sin(£;~a) 



We choose a normalized incident amplitude = 
l/y| sm(kfa)\. The net current primarily comes from 
transmission of electrons of energy between the elec- 
trodes' Fermi energies and is calculated as 0] 

/= tE / dET°{E)[f L {E)-f R {E)\ 

with the Fermi function f(E) = l/exp[(E - E F )/k B T] 
and the room temperature T — 300 K. For ferromagnetic 
electrodes the megnetoresistance is defined as the per- 
centage change of resistance between parallel and anti- 
parallel configurations R m = (iZ an « - Rparai) / R&nti = 

(Iparal I&nti) I Iparal • 

In metals, an electron of energy E may come from 
different energy bands. The corresponding effective pa- 
rameters e a m x an d Cx are t nen an average of the param- 
eters of these bands and are energy dependent. In lin- 
ear or quasi-equilibrium system, they are approximately 
the values near the Fermi energy. In a non-equilibrium 
system, if the difference of the Fermi levels between the 
two electrodes is comparable to the energy band width 
of the metals, the energy dependence should be taken 
into account. In the existing experiments the applied 
bias can be higher than 4 eV, which is larger than the 
width of the d bands where the Fermi level locates in 
many metals. In the case of Ferromagnetic Fe which ex- 
emplifies the electrode material here, approximately five 
bands can be identified from the density of states (DOS) 
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FIG. 3: (a) The I-V curve of a 30-basepair Poly(G)-Poly(C) 
DNA between two ferromagnetic Fe electrodes of parallel 
(solid curves) and anti-parallel (dotted) configurations. The 
thicker curves illustrate the total currents and the thinner 
ones the contribution from the two spin branches, (b) The 
magnetoresistance vesus the applied bias potential. In the in- 
set of (a), we show our theoretical fit (thick solid curve) to the 
experimental result (filled circles connected by dotted lines) 
in Ref. 0. 



near the Fermi energy of the bulk material • For the 
spin-up (majority) electrons, the five bands locate ap- 
proximately at 2.5, 0, —0.68, —3.4, and —7 eV above the 
Fermi energy with band width 6, 0.3, 0.6, 4.1, and 3.7 eV 
respectively. For the spin-down (minority) electrons, the 
energy bands are the same as above but shifted 2.58 eV 
to higher energy. Using Lorentzian broadening, we can 
mimic the bulk DOS and extract the parameters e m and 
t a m as shown in Fig. [2] At the Fermi energy, we get the 
hopping parameters 0.39 eV and 1.62 eV for spin-up and 
spin-down electrons respectively which coincide with the 
result obtained from the Fermi velocity (TJj ■ 

Just as in Ref. 0, we extract the parameters of the 
DNA chain by fitting the experimental result of Ref. 0] . 
By evaluating the energy-dependent par ameters and 
t a m from Platinum's band structure [l!j, we can fit the 
experimental result as shown in the inset of Fig.[3Ja). As 
a result, we find that the hopping parameter is td = 0.6 
eV, the equilibrium Fermi energy is 1.73 eV higher than 
the DNA HOMO on-site energy, the contact parameters 
are tdm = 0.019 eV for the right electrode and 0.013 eV 
for the left, 1/3 of the bias voltage drop at the right 
contact and 2/3 at the left. The above parameters are 
close to those obtained in Ref. except for a larger tdm 
in the present case. 

For ferromagnetic Fe electrodes as shown in Fig. EJa), 
we get two I-V curves corresponding to the parallel (solid 
curves) and antiparallcl (dotted curves) magnetization 
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FIG. 4: The magnetoresistance vesus the spin-flip parame- 
ter t a d in DNA when a bias of 4.8 eV is applied. The other 
parameters are the same as in Fig. [3] 



configurations for the left and right electrodes, using the 
same contact parameter td m = 0.02 eV for the two elec- 
trodes. In Fig. EHb), the magnetoresistance which de- 
scribe the percentage change of the resistance of the sys- 
tem when being switched from parallel to antiparallel 
configuration, is plotted. In contrast to the results where 
constant values of the parameters e a m and t a m at the Fermi 
energy were used, here we find a much smaller magne- 
toresistance until a strong bias is applied. This can be 
understood from the energy dependence of t m shown in 
Fig- Etc)- Instead of two parallel lines at 1.4 eV and 0.6 
eV, the curves of t m vs E for spin-up and spin-down elec- 
trons cross near the Fermi energy. This crossing makes 
the magnetoresistance disappear around the bias volt- 



age Vbias = 2 eV. The increase of magnetoresistance at 
higher bias voltage results from the increasing t m spread 
between spin-up and spin-down electrons in the range of 
2 eV around the Fermi energy. 

In Fig. we show how a spin-flip mechanism can affect 
the spin injection, assuming that the spin-flip can hap- 
pen only when an electron jump from one site to another. 
The magnetoresistance does not decay monotonically to 
zero. Instead, the magnetoresistance is enhanced when 
there is a very weak spin- flip coupling (t d < 1.9 meV) as 
a result of the quantum interference in the system. In 
the transmission spectrum, peaks are slightly split with 
the increase of t 8 d , indicating the mixing of the spin-up 
and spin-down states in the system due to the spin-flip 
coupling. We observe an increase in magnetoresistance 
from 20% at t s d = to 60% at t s d = 1.9 meV as displayed 
in the inset of Fig. Then the magnetoresistance de- 
creases smoothly until t d = 20 meV. Above that value, 
the magnetoresitance begins to oscillate when it decays 
to zero. 

In summary, we have investigated the quantum spin 
transport through a short DNA chain connected to ferro- 
magnetic electrodes. We have used a tight-binding model 
to describe the system where the parameters are ex- 
tracted from the experimental results and realistic energy 
bands of metals. We find that the energy band structure 
of the ferromagnetic electrodes significantly affects the 
resulting spin transport. In the presence of the spin-flip 
mechanism, enhancement and oscillation of magnetore- 
sistance due to mixing of spin states are also observed. 
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